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Abstract 

o . 

We present explicit formulas for the eigenvalues and eigenfunctions of the 
elliptic Calogero-Sutherland (eCS) model as formal power series to all orders in 
i-C ■ the nome of the elliptic functions, for arbitrary values of the (positive) coupling 

^p H , constant and particle number. Our solution gives explicit formulas for an elliptic 

| deformation of the Jack polynomials. 

B 

> 

1 Introduction 



The elliptic Calogero-Sutherland (eCS) system is a quantum mechanical model of iden- 
tical particles moving on a circle and interacting via a two-body potential given by the 
Weierstrass elliptic function p [HI ISul IQSj . It is denned by the 2nd order differential 
operator 

N r) 2 

H = -J2d^ + 7 £ v{ Xj -x h ) (i) 

j=l 3 l<j<k<N 

where N — 2, 3, . . . is the particle number, — n < Xj < n are coordinates on the circle, 

7 = 2A(A-1), A>0, (2) 

1 Contribution to the conference SPT 2004 in Cala Gonone (Sardinia, Italy) 
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is the coupling constant, and the two-body potential 



^4 sin [{r + i(3m)/2] 
which is essentially equal to the Weierstrass elliptic function p with periods 2tt and 

The eCS system is known to be integrable in the sense that there exist differential 
operators of the form 

N d n 
H n = > (— i) n — h lower order terms 

3 =l 3 

for all n = 1,2, . . . , N, which include the eCS Hamiltonian, H 2 = H, and which 
all mutually commute, [H n ,H m ] = for n,m = 1,2, ...,N |USj . Moreover, in the 
trigonometric limit (3 = oo where the two-body potential reduces to a trigonometric 
function, the explicit solution of this model was found by Sutherland more than 30 
year ago |Suj . In the two-particle case, N = 2, the eigenvalue equation of the eCS 
system is equivalent to the Lame equation studied extensively at the end of the 19th 
century; see |WWj for a review of the classical results. Recent work on the eCS model 
include jDH M lEFKl iFVll lFV2l IFHTl K]. 

In this paper we present a generalization of Sutherland's solution to the elliptic case 
without restrictions on parameters (see the Result in the final section). More specif- 
ically, we present explicit formulas for the eigenfunctions ^(x; 11 ) and corresponding 
eigenvalues £(n) of the eCS Hamiltonian, 

i7-0(x;n) = £(n)V>(x;n), x = (x x , . . . , x N ), (4) 

which are labeled by integer quantum numbers 

n = (n 1 ,...,n N ), n 1 >n 2 >...>n Nl (5) 

and which are of the following form, 

■0(x; n) = $(x; n)\P(x) (6) 

where 

¥(x)= J] e( Xj -x k ) x (7) 

l<j<k<N 



2 To be precise: V(z) = p(z) + c with c = 1/12 - (1/2) J2m=i sinh~ 2 [(/3m)/2]. 
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with 

oo 

e{z) = sin(z/2) JJ(1 - 2q 2n cos(z) + g 4n ) , q = e"^ 2 (8) 

n=l 

essentially the Jacobi Theta function $i. 3 The $ are symmetric functions of the vari- 
ables zj = e 1Xj , and they are in one-to-one correspondence with the plane waves 

JV 

S(x;n) = ]T JJeNi-i (9) 

Q£S N j=l 

which provide a complete set of eigenfunctions in the non-interacting case 7 = 0. It is 
important to note that our solution in the trigonometric limit reduces to Sutherland's: 
the eigenvalues £ (n) become equal to the well-known expressions 

N 2 

£ (n) = J2(n 3 + X[l(N + l)-j]) , (10) 

3=1 

and the functions $(x;n) reduce to the Jack polynomials playing a prominent role 
also in various other contexts in mathematics; see jMcDl lSt]. We also note that, for 
(3 < oo, the functions $(x; n) are no longer polynomials, and also the eigenvalues £(n) 
become much more complicated. In particular, \1> is not the ground state of the eCS 
Hamiltonian for (3 < oo. Correspondingly, our solution is by infinite series which, 
at this point, are only formal: we leave open the important but difficult questions 
of convergence and resonances (as discussed in more detail below). We only mention 
the results in Ref. |KTj which suggest that our series solutions have a finite radius 
of convergence in the nome q of the elliptic functions, and our results suggest that 
there exists a resummation so that the resonances disappear |L4j . Moreover, in the 
trigonometric limit resonances do not appear |L2j . and all our infinite series collapse 
to finite ones. 

The present paper is based on our previous results in |CI4 ILH IL21 IL31 IL4j . Our 
starting point is a Theorem obtained by quantum field theory techniques in |L2j and 
proven by direct computations in |L4j . This theorem suggests to write the eigenfunc- 
tions $(x;n) as series of particular symmetric functions F(x; n) which are given by 
the following explicit formulas, 

cv \ Till ^3 t n il ni<j<fc<jv ®(£j/£fc) A 
F ( x ' n ) = li|£ 2rt*i' tin -^fc^x ( U ) 



n;: fc=1 0(e-ve 



3 To be precise: 0{z) = Mz/^)/[^Q 1/4 Un=i( l - 1 2n )\- 
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where 

oo 

0(0 = (i - II K 1 - ^0(1 - 9*70], (12) 

m=l 

and the integration contours are nested circles in the complex plane enclosing the unit 
circle, 

c: ; :( <•'<••'• -7r<yj<TT, 0<e<f3/N. (13) 

Note that the F(x; n) are symmetric functions of the variables Zj = e 1Xj , and they 
can be expanded as Laurent series. A simple but important consequence of the above 
mentioned theorem is the following. 



Corollary: Let 

$(x;n) = ^a(/x;n)F(x;n + /x) (14) 

m 

where the sum is over all 

fx= 22 Vjk^jk with n jk G Z and (E jk ) £ = 5 je - 5 ke (15) 

l<j<k<N 

for £=1,2,..., N, and 

a(M;0) = 5(/x,0)+O( 7 ). (16) 

Then n) defined in Eqs. (0) and |7|) is an eigenfunction of the eCS Hamiltonian 
with corresponding eigenvalue £(n) provided that 

[£ (n + /£) - £{n)]a(fj/,n) = S u a(/j, - iyE jk ; n) (17) 

j<k uei* 

with ^o(n) in Eq. MU\) and 

v isq 2l/ 

S = 0, S v = =- and S_ v = — for v > 0. (18) 

1 — q iv 1 — q AV 



Note that S v = [1 — 5(z/, 0)]z//(l — g 2iy ) for all integer v, but we prefer our somewhat less 
elegant definition which makes manifest that, in the trigonometric limit q = 0, S v = 
for v < 0. This implies that, for q = 0, Eq. ()17j) has triangular structure (in a natural 
sense explained in Ref. |L2J), which implies that the eigenvalue £(n) is equal to £o( n )- 
Moreover, one obtains a simple recursion relation to compute all the non-zero a(fi; n) 
in in Eq. (|14|) . It is interesting to note that this solution algorithm for q = is different 
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from Sutherland's, even though is yields the same solution |L2j . We did not realize in 
|L2j that it is possible to obtain an explicit solution of the recursion relations for the 
coefficients a (/a; n) as follows, 

oo oo 
s=0 ji<ki js<k s v\,...,v a =l 

x— Ffi — tt^f ^ v— forg = 0, (19) 

with 

6(/x;n) = £ (n + ^) -£ (n), (20) 

where all sums are in fact finite due to the Kronecker delta and certain properties 
of the functions F(x; n) discussed in [L2J (Eq. (fT9*|l is a simple special case of our 
general result presented in the last section). This together with the results in Ref. 
|L2j provides explicit formulas for the Jack polynomials. It is interesting to note that 
the recursions relations which one gets in Sutherland's algorithm are more complicated 
and, to our knowledge, have not been solved explicitly for general N (we are only aware 
of similarly explicit previous results for N = 3 jPRZj ). It is also worth noting that, due 
to translation invariance, the dependence of the eigenfunctions on the center-of-mass 
coordinate x\ + • • ■ + £jv is trivial, and therefore the functions $(x; n) are of the form 
e mjv(a;i+-a!jv)$^ x . w here hj = nj —tin obeys h± > rii > . . . njv-i > = 0. However, 
it seems that this is not manifest in our explicit formula for these functions, and we 
therefore seem to get an infinite number of different representations (labeled by njv) 
for each distinct Jack polynomial (this is true not only for q = but also in the elliptic 
case). 

In the rest of the paper we discuss how to explicitly compute the coefficients a(/x; n) 
and eigenvalues £(n) for the general elliptic case from the Corollary above. The strategy 
is to find a "good" expansion parameter, allowing to solve Eq. (|17|) recursively. The 
obvious parameter is the squared nome q 2 of the elliptic function LI . While this 
provides a possible solution algorithm, we were able to obtain the explicit solution 
only up to order (q 2 ) 7 for N = 2 with the help of MAPLE in this way |L4j . The 
formulas obtained are rather complicated and suggest that it is hopeless to find explicit 
expressions to all order in q 2 . However, this result motivates a more efficient solution 
strategy by expanding in the coupling parameter 7. As we show, the resulting solution 
algorithm is indeed simpler, and while we obtained the explicit solution up to order 
7 6 without the help of MAPLE, it again seems hopeless to obtain explicit formulas 
to all orders in 7 in that way. However, this result gives a better understanding of 
the structure of the solutions, and it led us to a method of solution to all order. The 
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key to this was to introduce a parameter r] "by hand" which efficiently organizes the 
complexity of the solution, and this allows us to obtain the solution as a power series 
in 77 to all orders (without the help of MAPLE). From this our explicit formulas for 
£(n) and a(/x; n) of Eq. (fTTjl are obtained by setting 77 = 1. 

In the next three sections we outline the explicit solutions obtained by expanding 
in q 2 , 7 and 77, respectively. For simplicity in notation we restrict this discussion to the 
simplest non-trivial case N = 2. The generalization to arbitrary iV is straightforward, 
but we only give the result in the final section. We plan to include a more detailed 
derivation of this in a future revision of Ref . |L4j . 

2 Expanding in q 2 

For simplicity we restrict ourselves to iV = 2. In this case [i = /iE 12 , and (fTTj) simplifies 
to ~ 

[£ (n + fiE 12 ) - £ (n)]a(//;n) = 7 y~] S v a{n - v\ n) (21) 

where a(/z; n) is short for a(/xEi 2 ; n). To simplify notation we suppress the dependence 
on n in the following. We make the ansatz 

00 00 

a(ii) = ^(aO?*, £ = Y. S M 2e > ( 22 ) 
1=0 1=0 

and with 

ot/(0) = 6(p, 0) and a e (fi) = for \i < -i (23) 
we obtain by simple computations (expanding the S v in geometric series etc.) 

t i 
b(fjL)a£t(fi) = ^2 £m(xe-m(fJ>) + 7 X] va ^ - v ) 

m=l v=\ 

I t/v 

+7 X] v \ a ^m{^ -v) + ae-um(n + v)\ (24) 

v=l m=l 

where b(fi) = £ (n + /ttE 12 ) - £o(n). Using £ (n) = (m + A/2) 2 + (n 2 - A/2) 2 we get 

b(ji)=2fi(P + fi), P = n 1 -n 2 + X. (25) 

It is important to note that Eq. (|24|) has triangular structure: for each i and /x 7^ it 
determines a^(/i) as a finite sum of terms involving only a^(/i') and £^ with £' < £ and 
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ottifjf) with [t! < /i, and the equation for \i = allows to determine the Ei recursively. 
We also note that the conditions in Eq. ()23|) is essential for getting a simple recursion 
procedure. It is straightforward to implement the recursive computation of the ae(n) 
and £i in a symbolic computing software like MAPLE or MATHEMATICA. We only 
give the results for the eigenvalues which we obtained using MAPLE, 
1 



£1 



P 2 -i 



7 
1 



Sa 



(P 2 -4)(P 2 -1) 
1 

(P 2 -9)(P 2 -1) 



1 

(P 2 -16)(P 2 -1) 



6(P 2 - 2) 7 2 - 6 7 3 + 



(5P 2 +7) 4 
4(P 2 -1) 2 ' 



12(P 2 - 3)7 2 -487 3 + 



4(15P 4 -37P' 



4(7P 2 



-112 7 



(P 2 -4)(P 2 -1) 

(7P 4 -74P 2 +112) 
4 (P 2 -4) 

3(365P 10 -6662P 8 - 



(P 2 -4)(P 2 - 

(9P 4 +58P 2 +29) 6 
2(P 2 -4)(P 2 -1) 4 7 



l) 2 ' 



(26a) 
(26b) 

(26c) 



7 2 - I8O7 3 

f42249P 6 -115640P 4 +119816P 2 - 



2(P 2 -9)(P 2 -2) 2 (P 2 -1) 2 

3(259P 8 -3358P 6 +11415P 4 -4252P 2 -25664) 
(P 2 -9)(P 2 -4) 3 (P 2 -1) 2 



7 J 



2151P 10 -18127P 8 -10529P 6 +293115P 4 -501962P 2 - 



715P 8 



4(P 2 -9)(P 2 -4) 3 (P 2 -1) 4 
-481P 6 -43203P 4 +94061P 2 +104428 



■18528) 4 



79832 ^,6 



(26d) 



4(P 2 -9)(P 2 -4) 3 (P 2 -1) 4 



7' 



I 1469P 10 +9144P 8 -140354P 6 +64228P 4 +827565P 2 +274748 8 
i " 64(P 2 -9)(P 2 -4) 3 (P 2 -1) 6 ' 

We computed £^ up to order 1 = 1 but for obvious reasons do not write down our 
full result. This result was previously obtained in [FGP] up to t = 2 using a different 
method, and we convinced ourselves that our results agree. 



Remark: In the computations discussed above we implicitly assumed that 6(/i) is 
always different from zero, which is only the case if A, and thus P, is not an integer. If 
6(/i) vanishes we have a resonance, and while it is possible to generalize the algorithm 
to allow for resonances |L4| we will not discuss this here for simplicity: Throughout 
this paper we ignore resonances. (As discussed in |L4j . we believe that resonances are 
not a serious problem.) .) 



We observe that the Ee become more and more complicated with increasing £, and 
it seems that they are polynomials in the coupling constant 7 as follows, 

s=2 
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(s) 

where E\ become more complicated with increasing s. It is interesting to note that the 

(s) 

formulas for the coefficients E\ become somewhat simpler when expanded in partial 

(2) 

fractions. In particular, we found by inspection that all E\ we computed can be 
written in the following simple form, 



V M l —\ (2f 



(the sum is over all integer divisors k of £) . We checked this formula up to £ = 7, and 
by assuming it to be true for all i we obtain by a simple computation 



OO OO / 

£=0 k=l ^ 



1 ^ ^ (») 



k P + kj 2(1 - q 



(this conjecture will be proven in the next section). This formula suggests that it is 
possible to obtain explicit expressions to all orders in q 2 if one expands in the coupling 
parameter 7. We now present an alternative and simpler solution algorithm motivated 
by this observation. 



3 Expanding in 7 

We now make the ansatz 

i(A*) = 5> ( *V, £ = S + J2£ [s h s , (30) 



00 00 

'»(, 

s=0 s=l 



and by simple computations we obtain from Eq. (JTTj) 

s 

b{fi)a {s) (n) - £ {s ~ r) u ir) (m) = S » a{S ~ 1] (V ~ v ) ( 31 ) 
which we can solve with the following ansatz 



a 



(0) 



(n)=6(ji,0), a (s) (0) = for s = 1,2,.... (32) 



This yields, in particular, 

S {s) = -^^« (s - 1} (-z/). (33) 
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The recursion relations in Eqs. (j3*Tj) are much simpler than the ones in Eq. f!24|) . By 
inspection we find that the following ansatz is consistent, 



a {s) (fi)= ^2 S v ,---S n f a (y 1} ...,u a ;ii), (34) 
ui,...,u a ez 

and by straightforward computations, 

fiWlV = u \ (35a) 

fzivuvXM = 77 : yrr^ 35b 

M^I^^M) = 77 777 T 777"^ 77 777 ^ 2 (35c) 



/ 4 (z/i,I/ 2 ,^3,^4;/i) 



5(/i, Z/i + f 2 + f 3 + ^4) 



6(z/i + z/ 2 + z/ 3 + ^4)^(^1 + ^2 + vz)b(vi + ^ 2 )fe(^i) 



b{y 2 )b{y x fb{y x + i/ 4 ) 6(z/ 4 )&fa + ^) 2 &K 

<K^2 + ^3 + ^4,0)5(/i, Z/i) 



6(^3 + V4)b(v 3 )b(vi) 



(35d) 



etc. (it is useful to note that the variables Uj can be permuted in each term). Using 
Eq. (j22D we get 

£ {s) = - s v.---S Ul fs-i(yi,---,Vs-i\-Vs), (36) 

in particular, £^ = and 

5 (2)_ y-c o ^1 + ^2,0) _ y^ 55 (J_ + _L_\ (37) 

which proves Eq. ([29)1 . With the formulas given above one can easily write down 
similarly explicit formulas for 8^ for s = 2, 3, 4, 5. 

Computing the functions f s up to s = 6 we observed the following simple patterns: 
it seems that the building blocks for the solution are the following expressions, 

Gf = £ S,....^ f; ^ 1 + ... + fca _ 1 -£) (38) 
^,...,i/ s ez fei,...,fc s _i=o YYr=\°\l^t=\ v e-) 
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With that we can write the coefficients of the eigenvalues in a simple manner as follows, 





— — Ltq 


(39a) 


£(3) 


= -G?> 


(39b) 


£(4) 


r (A) r (2) r {2) 


(39c) 


£(5) 


— — Ltq + Ltq Ct x + Cr Q Lr x 


(39d) 




^(6) , ^(4)^(2) ^(2)^(4) 






,^(3)^(3) ^(2)^(2)12 rr ,(2) l2r( (2) 


(39e) 



This suggests that 

oo 

£ (s) = £(-l) P (- • O^i + • • • «n + • " • + *r, r - 1)GW • • ■ g£° (40) 

r=l 

where '(■••)' are possible combinatorial factors which might appear at higher order 
but, up to s = 6, all are equal to 1. 

We find that, for s > 6, nontrivial combinatorial factors in the formula above 
appear, and this destroys the hope that we can find a closed formula for £ ^ for 
arbitrary s in this way. Still, Eq. (|4()|) suggest that we can write £ in a simple manner 
using the following quantities, 

oo 

G fe = J>£V, (41) 

s=2 

and the formulas above suggest, 

£ = £ - G + God - G [(?i] 2 - [G ] 2 G 2 + ... (42) 

where the dots are higher order terms. We will prove and extend this formula in the 
next section: We will obtain an explicit formula of the following kind, 

£ = £ + J2^ £>= Yl (■■■)G n ---G rn (43) 

n=l r\,...,r n 

with certain combinatorial factor '(•••)' which we will compute explicitly. It is interest- 
ing to note that £ n is of order q 2n (since all G s are of order q 2 ), and thus this formula 
allows to deduce in a simple manner the series expansion in q 2 . However, the term 
£ n contributes to all the (g 2 ) m -terms m > n, which explains why our expansion in q 2 
yielded so complicated expressions. We stress that, since £ n = 0(q 2n ), Eq. (|50J) still is 
an expansion the q 2 and has, as we believe, a finite radius of convergence. 
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4 Expanding in 77 



The results in the previous section led us to a more powerful solution strategy which 
we now explain. Defining an operator § as follows, 

S/M := -i/) (44) 



we can write Eq. (|17|) as 

where 

Making the ansatz 



[&0) - £]a(/i) = 7§«(/i) (45) 
5 = £ - £ . (46) 



a(/x) = 5(^0) + ^a s (/i) 7 s (47) 

s=l 

we get [b(/i) — £]a s (fi) = §a s _i(yu), and thus 

oo 

«(A*) = E7*([60')-^]" 1 S) 5 (5(M,0) 

oo 1 



s=0 

1 1 



=0 ^.....^.ez HfJ-) ~ £ ) ~ v a ) - £} [b(p -v a - u a -!) - £} 

1 



[&(// - i/ s - £} 

oo 



~-6(fJ, -v a - U s -! l/i,0) 



£7* £ ■V---S-. r J^ lVr) z J V (48) 



Setting /i = in Eq. (|4*HJ) gives 5 = — 7§a(0), which implies 

£-_yV +1 V 9 • ••<? <HEr=l^,Q) r4 x 

It is easy to see that the s = term here vanishes, and by a shifting the summation 
variable we obtain the following equation determining £, 

£ = -G{£) (50) 
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where 



G (o ■= E^ s E s ^--- s ' 



"n::;[KELi^)-e] 



s=2 fi,...,i/ s GZ 

We now observe that the functions G(£) has a Taylor expansion as follows, 

oo 

g(o = E 



(51) 



(52) 



fc=0 



with 



1 



E^ E ^-s 



^(e:=i^,o) 



s=2 i/ lv ..,i/ s 



E^ s E Il s » E *(*>£U*r) 

s=l i^,...,i/ s £Zr=l fci,...,fc s =0 



^(e: =1 ^ ; o) 



n^iP»(E5=i^)-e] 



i+k s 



(53) 



Note that these are exactly the quantities which we found empirically in the last Sec- 
tion. To solve Eq. (}50|) efficiently we replace it by 



£ 



-Q 



(54) 



k=0 



where we introduce a parameter rj serving as useful book keeping device (we set r\ — 1 
at the end of the computation). It is straightforward to compute the Taylor series of 
£ recursively: With the Ansatz 



£ = J2^v k 



(55) 



k=l 



we get by simple computations 

£\ — —Go 
£2 = GoG\ 

£?. = — [Gq] 2 G2 — Gq[Gi} 2 



-3 
£4 



?>{Gq\ 2 G\G2 + Got^i] 3 + [Go] 3 G3 



£ 5 = -4[Go] 3 G 1 G 3 -6[Go] 2 [G 1 ] 2 G 2 -2[Go] 3 [G 2 ] 2 -[Go] 3 G4-Go[G 1 ] 4 (56) 
etc. This suggests that 



fc ,fci,---,fc„_i=0 



^ = ("I)" E - !> T.U ^ ( ko, U . • , in- I 



n-1 



3=0 



(57) 
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for all n = 1, 2, . . . (we checked that using MAPLE up to n = 15). This result is, in 
fact, a classical theorem due to Lagrange: 4 the equation determining £ is of the form 
8 = —r)G(E), and thus Lagrange's theorem as stated in jWWj. Paragraph 7.32, implies 



-v) 



n=l 



III 



n An— 1 

dy n-i ^ 



(5? 



equivalent to Eq. (|57j). 



5 Conclusions 



It is possible to generalize the results of the previous two sections to arbitrary particle 
numbers N. We intend to give the details in a future revision of Ref. |L4j and quote 
here only the result. 



Result: The eigenvalues of the elliptic eCS model are given by 

oo oo 

£(n)=£b(n)+X;(-l) B £ 

n=l feo,---,A; n _i=0 
/ 1 \ n ~ l 

xd (n — 



'An - 1, EJlJiWn, El:! 1 *i) L " . ) ^ J II[Gi(n)]^ (59) 
wif/i £o( n ) = Y^f=i (jtj + (-^ + 1) — j]j ! n \ >U2> ■ ■ ■ un integers, and 

oo oo 

G*(n) = £f£---£ E <v-&. £ 

s=2 ii<fci j s <fc s yi,...,v s ez h,—,£s=o 

xS(k,T, S r x l r ) ^EU^^O) 

n r =i KEw^ E i<*«i n ) * 

w/jere 6(/x; n) = £o(n + /x) — £o(n) ; = [1 — <5(z/, 0)]i//(l — g 2l/ ) and E jfc m #g. (T73)) . 
T/ie corresponding eigenf unctions are given by Eqs. (0J), (0) and Jj^| )-(f73 |) with the 
coefficients 



s=i ii<fci is<fc s v!,...,v s ez 



4 I thank S.G. Rajeev and G. Lindblad for helpful discussions on this. 
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Kjb Er=l ^EjyfeJ 

n: =1 [KELi^n)-£(n)+£ (n)]- 



(61) 



There is an even more explicit formula for the coefficients a{n) which we plan to 
give elsewhere. 

We do not label this as theorem since we did not check the details of our formulas 
here carefully enough to be sure that there are no typos and/or (minor) mistakes. 5 
The purpose of this paper was to make available the explicit solution of the eCS model 
which we announced in two recent meetings, and to describe the last part of a somewhat 
lengthy journey leading us to this result. We feel that we are not at the end of this 
journey yet: quite some work remains to be done to understand this result. 

Anyway, it seems fair to say that the eCS model is an exactly solved model now. 
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